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A  RESUME  OF  DETERMINISTC  TIME-VARYING  LINEAR  SYSTEM  THEORY 
WITH  APPLICATION  TO  ACTIVE  SONAR  SIGNAL  PROCESSING  PROBLEMS 


by 

Lewis  Meier 


ABSTRACT 


Brief  summaries  are  given  of  time-varying  linear  system  theory  in  both  the 
convention  frequency-shift  version  and  the  novel  frequency -shift  version. 
The  latter  version  of  the  theory  is  utilized  to  obtain  the  matched  filter 
output  for  the  bistatic  echo  from  a  moving }  turning  line  target  for  a 
linear  FM  pulse  in  a  non-dispersive  medium. 
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INTRODUCTION 


The  major  concern  of  this  report  is  with  a  sonar  return  processed  in  the 
standard  manner  with  a  beam  former  (not  considered  herein)  and  a  bank  of 
matched  filters,  where  the  modification  of  the  transmitted  signal  by  the 
medium  and  the  target  can  be  represented  by  a  time-varying  linear  system. 
Two  transmitted  signals  are  of  interest:  the  CW  pulse,  which  has  unity 
time/bandwidth  product  and  the  FM  pulse,  which  has  a  time/bandwidth  product 
much  greater  than  unity.  Correspondingly  this  report  summarizes  two 
versions  of  deterministic  time-varying  linear  system  theory:  the  standard 
frequency-shift  version,  applicable  to  analysis  of  detection  and  measure¬ 
ment  extraction  using  CW  pulses  and  the  more  novel  frequency  stretch 
(doppler)  version,  applicable  to  analysis  of  detection  and  measurement 
extraction  using  FM  pulses. 

The  approach  is  to  characterize  the  various  parts  of  the  signal  path 
between  a  sonar  transmitter  and  a  sonar  receiver  as  linear  time-varying 
systems.  The  theory  then  provides  means  to  characterize  the  signal  at 
various  point  along  the  path,  to  characterize  the  parts  of  the  path 
between  these  points  and  to  relate  these  characterizations.  Specifically 
the  signals  are  characterized  by  cross-ambiguity  functions  (called  the 
ambiguity  function  for  the  transmitted  signal),  which  tell  what  the 
response  to  the  signals  of  a  bank  of  matched  filters  would  be  as  a  function 
of  time  shift  and  either  frequency  shift  or  frequency  stretch;  the  parts 
of  the  path  are  characterized  by  spreading  functions,  which  tell  how  the 
parts  spread  the  signal  in  time  shift  and  frequency  shift  or  frequency 
stretch  and  the  relationship  is  a  modified  double  convolution  in  time 
shift  and  frequency  shift  or  frequency  stretch. 

Rihaczek  [l]  made  an  extensive  study  of  the  use  of  matched  filter  radars 
and  since  he  was  primarily  interested  in  point  targets  and  a  perfect 
medium,  he  investigated  in  great  detail  ambiguity  functions  of  the 
transmitted  signal.  In  fact  his  book  could  be  titled  "All  You  Wanted  to 
Know  about  Ambiguity  Functions  but  were  Afraid  to  Ask".  The  use  of  the 
standard  theory  was  applied  to  sonar  problems  by  S0strand  [2],  whose 
development  is  followed  here  with  minor  modifications.  Kelley  and  Wishner 
[3]  studied  generalization  of  the  classical  ambigu'ty  function  including 
the  doppler  ambiguity  function  studied  herein.  They  were  also  interest  in 
the  situation  in  which  acceleration  of  the  target  caus  d  a  significant 
change  in  velocity  during  the  duration  of  the  transmitted  pulse.  While 
acceleration  can  be  very  important  with  targets  such  as  m;  riles,  it  seldom 
is  so  in  the  sonar  situation  and,  hence,  is  neglected  herein.  The  complete 
development  of  the  frequency-stretch  theory  presented  here  is  believed  to 
be  novel . 

Weston  [ 4 ;,  Kramer  [5]  and  Russo  and  Bartberger  [6]  have  all  derived  the 
doppler  ambiguity  function  of  a  linear  FM  pulse  in  terms  of  Fresi  el 
Integrals  with  differing  limits  of  integration.  Russo  and  Bartberger  give 
the  exact  limits,  which  are  rather  complicated,  while  Weston  uses  a  gross 
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approximation  to  obtain  the  general  behaviour  of  the  doppler  ambiguity 
function  not  only  for  linear  FM  pulses,  but  curved  FM  pulses  as  well. 
Weston's  paper  is  highly  recommended  for  its  physical  insight  into  the 
nature  of  ambiguity  functions  of  FM  pulses.  Kramer  (and  T.  Kooij  in 
unpublished  work  of  1964)  uses  limits  of  integration  that  are  valid 
approximation  for  time  shifts  small  compared  to  the  pulse  lengths;  herein 
the  next  level  of  approximation  to  the  limits,  which  allows  large  time 
shifts  but  ignores  changes  in  length  of  the  signal,  is  derived.  This 
approximation  is  valid  for  normal  sonar  situations  when  the  time-bandwidth 
product  of  the  FM  signal  is  less  than  750.  Furthermore  it  is  shown  that 
when  in  addition  the  time-bandwidth  product  exheeds  50,  the  Fourier 
transform  with  respect  to  time  shift  of  the  doppler  ambiguity  function  of 
a  linear  FM  pulse  takes  a  particularly  simple  form. 

A  good  and  useful  example  of  the  application  of  the  time-shift  theory 
presented  herein  is  determination  of  the  matched  filter  output  for  the 
bistatic  echo  from  a  moving,  turning  line  target  for  a  linear  FM  pulse, 
when  all  medium  effects  other  than  propagation  delays  are  ignored.  The 
doppler  spreading  function  of  such  a  target  may  be  readily  approximated 
in  a  simple  form  and  then  used  to  obtain  the  matched  filter  output  as  a 
convolution  between  an  impulse  response  function  representing  the  target 
and  a  window  function  representing  the  FM  pulse.  This  result,  derived 
herein,  is  also  believed  to  be  novel  and  provides  the  basis  in  a  sequel 
for  algorithms  for  measuring  target  properties. 

Two  topics  I  will  not  be  concerned  with  in  this  report  are  propagation 
loss  and  background.  The  sonar  equation  and  its  extension  may  be  utili¬ 
zed  to  determine  level  of  the  echo  with  respect  to  background  at  the 
matched-filter-bank  input;  I  am  concerned  with  how  that  echo  is  spread 
in  time  and  frequency  or  Doppler  at  the  matched-filter-bank  output. 

Both  pieces  of  information  are  required  for  investigating  detection  and 
measurement  extraction.  The  theory  presented  herein  and  in  a  sequel  on 
stochastic  systems  have  direct  application  to  the  analysis  of  reverbera¬ 
tion,  but  I  will  not  pursue  that  topic  either  herein  or  in  the  sequel. 

The  main  results  are  contained  in  the  body  of  the  report  with  the 
detailed  derivations  relegated  to  appendices.  The  body  of  the  report  is 
divided  into  two  sections:  the  first  is  concerned  with  development  of  the 
theory  while  the  second  is  concerned  with  its  application  to  finding  the 
matched  filter  output  of  the  bi static  echo  from  a  moving  turning  target 
for  a  linear  FM  transmission  in  a  non-dispersive  medium. 
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1  THE  THEORY 

Figure  1  is  a  simplified  block  diagram  of  a  single  beam  of  an  active  sonar 
system.  The  matched  filters  bank  consists  of  a  set  of  filters  matched  to 
the  transmitted  signal  x(t)  shifted  in  frequency  by  different  amounts. 

In  practice,  for  a  CW  pulse  the  matched  filter  bank  is  realized  by  taking  a 
time  sequence  of  Fourier  transforms  of  segments  of  echo  whose  length  approx¬ 
imately  equals  that  of  the  pulse,  while  for  an  FM  pulse  only  a  single 
matched  filter  matched  to  the  transmitted  signal  is  used  because  of  the 
large  doppler  tolerance  of  this  signal.  This  form  of  processing  is  predi¬ 
cated  on  the  assumption  that  the  doppler  effect  is  a  frequency  shift,  when 
in  actual  fact  it  is  a  frequency  stretch;  therefore,  the  matched  filter 
bank  should  consist  of  a  set  of  filters  matched  to  x(t)  stretched  in 
frequency  by  different  amounts.  Fortunately  for  an  CW  pulse  the  frequency 
stretch  is  adequately  approximated  by  a  frequency  shift  while  for  an  FM 
pulse  use  may  be  made  of  the  equivalence  between  zero  frequency  shift  and 
unity  frequency  stretch. 


FIG.  1  SIMPLIFIED  BLOCK  DIAGRAM  OF  A  SINGLE  BEAM 
OF  AN  ACTIVE  SONAR  SYSTEM 

1 .1  The  cross-ambiguity  function 

Since  we  are  trying  to  analyze  the  matched-filter-bank  ouput,  it  is 
reasonable  to  characterize  the  signal  y(t)  at  various  points  between 
transmission  and  reception  by  what  output  it  would  yield  if  used  as  input 
for  the  matched-filter  bank.  This  function  of  the  time  shift  t  and  the 
frequency  shift  v  or  the  frequency  stretch  a  is  formally  known  as  the 

cross-ambiguity  function  <j>.  or  <{>?.  Throughout  the  report  the  super¬ 
script  D  —  for  doppler  —  will  be  used  to  denote  the  frequency  stretch 
version  of  variables.  In  Fig.  1,  the  cross-ambiguity  function  <f>  or 

D  3 

<J>  of  the  echo  is  the  matched-filter  bank  output,  while  the  ambiguity 

3  D  D 

function  of  the  input  y  =  <Po  or  y  =  <j>Q  (where  the  cross  prefix  is 

dropped  for  obvious  reasons)  tells  how  the  transmitted  signal  is  spread  in 

x  and  v  or  a  .  The  cross-ambiguity  functions  <J>.  and  4>9  are  defined 

by  1  1 

00 

4>i(x,y)  =  j  e-2irjv(t-x)  x*(t_T)  y.(t)  dt  ,  (Eq.  la) 

-00  oo 

<t>i  (x,a)  =  a  j  x*[a(t-x)j  y^t)  dt.  (Eq.  lb) 

o 

where  *  denotes  complex  conjugate. 
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1.2  The  spreading  function 

The  effect  of  the  target  on  the  signal  is  to  shift  it  in  time  and  doppler 
due  to  its  range  and  range  rate  and  to  spread  it  in  time  and  doppler  due 
to  its  apparent  length  and  apparent  turning  rate.  It  is  of  course  these 
effects  that  can  be  used  to  measure  these  properties  of  a  target  as 
Plaisant  has  so  ably  investigated  [7,8].  On  the  other  hand,  the  effect  of 
the  medium  on  the  signal  is  to  spread  it  in  time  and  frequency  due  to 
multipaths  and  motion  in  the  medium.  Since  a  time-varying  linear  system 
such  as  a  target  or  the  medium  modifies  its  input  yj(t)  in  time  and 

frequency  or  doppler;  it  is  not  unreasonable  to  suppose  that  its  output 
y^(t)  be  a  weighted  sum  of  the  input  shifted  in  time  and  either  shifted  or 

stretched  in  frequency: 

y0(t)  =  _/  /j^^'^yjU-i)  Sqj (t ,v)  dvdx  ,  lEq.  2a) 

=  f  J  yjalt-x)]  Sj0(i,a)  dvdx  .  (Eq.  2b) 

-ao  o 

Equations  2a  and  2b  define  the  spreading  functions  3qj(x,v)  and 
Sqj(x,cx)  of  the  system  producing  yQ( t )  from  yj(t). 

It  is  more  standard  to  describe  a  time-varying  linear  system  by  its 
impulse  response  hjQ(t,x),  which  is  defined  as  the  system  output  at  time 

t  due  to  an  impulsive  input  at  time  t-x.  Use  of  an  impulsive  input 
yj  =  6(t-x)  in  Eqs.  2a  or  2b  yields: 


h 


ioU.t) 


=  ;  e2irJv(t-x) 


Sjq(x.v)  dv 


J 


=  /  I 

o  a 


<-D  <  (a-1 )  t  +  x 
^IOl  a 


,  a)  da  . 


(Eq.  3a) 


(Eq.  3b) 


From  the  similarity  of  Eq.  3a  to  a  Fourier  transform  it  is  clear  that  it 
can  be  inverted  to  yield 

SI0(x,v)  =  7  hIQ(t,x)  dt.  (Eq.  4) 

By  its  very  linearity  every  time-varying  linear  system  has  an  impulse 
response;  therefore  it  has  also  an  ordinary  spreading  function.  On  the 
other  hand,  it  is  apparent  from  the  impossibility  of  inverting  (Eq.  3b) 
that  not  every  time-varyinq  linear  system  has  a  Doppler  spreading  function. 
Nonetheless  some  very  important  systems  do  have  Doppler  spreading  functions. 
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1 . 3  The  fundamental  relationship 

Of  key  importance  to  the  analysis  is  the  relationship  between  the  cross- 
ambiguity  functions  <j>j  or  <j>j  and  4>q  or  <t>jj  of  the  input  and  output 
of  a  linear  time-varying  system  as  characterized  by  its  spreading  function 
SjQ  or  Sjq  ;  this  relationship  is  a  modified  convolution  in  t  and  v 
or  a  as  illustrated  in  Fig.  2.  By  applying  this  result  recursively 
starting  with  <j>  or  <f>^  (the  ambiguity  function  of  the  transmitted 

0  r\  0  r>  *  p 

signal),  <j>i  or  4>i  ,  <p2  or  <p2  and  finally  <j>3  or  <j>3  (the  cross- 
ambiyuity  function  of  the  echo)  may  be  computed.  Alternatively  the 
relationship  between  the  overall  spreading  function  S^G  or  sjjc  of  the 

two  systems  with  spreading  functions  or  S^B  and  SgC  or  SgG, 

which  is  the  same  modified  convolution  in  t  and  v  or  a  as  is 
illustrated  in  Fig.  3,  may  be  applied  recursively  to  find  the  spreading 
function  between  transmitted  signal  and  echo  and  then  the  original 
relationship  between  transmitted  signal  and  echo  used  to  find  <j>3 

or  $?. 


■  *1  ► 

SIO 

c 

^0  =  SIOT,V^I 

. . .  .  ^ 

or  y>P 

I 

D 

or  Sjo 

D  _D  A 

^0  =  Si0T,a^i 

FIG.  2  RELATIONSHIP  BETWEEN  THE  INPUT  AND  INPUT  CROSS-AMBIGUITY 
FUNCTIONS  AND  THE  SPREADING  FUNCTION  OF  A  TIME-VARYING 
LINEAR  SYSTEM 


FIG.  3  THE  SPREADING  FUNCTION  OT  TWO  SYSTEMS  IN  TAMDEM 
IN  TERMS  OF  THEIR  INDIVIDUAL  SPREADING  FUNCTIONS 
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I  I 

The  modified  double  convolution  operators  (x)  and  (x) 
defined  by  t,v  x,a 


required  are 


b  (x) '  a(i,v)  =  /  /  e^v  T  ^  a(x-x',v-v')  b(T',v')  dv1  dx'  , 

-OC  —  CO 

T,V 

(Eq.  5a) 


bD  (x)  '  aD(x,a)  =  /  /  a^[a‘  (x-x 1 ) ,  b^(x',a')  da'  dx1  . 


x  ,a 


"00  0 


(Eq.  5b) 


By  rewritting  the  definitions  (Eqs.  la,  lb,  2a  and  ^b)  in  terms  of  these 
modified  convolutions,  it  may  be  shown  that  the  relationships  depicted 
in  Figs.  2  and  3  hold. 


1.4  Some  approximate  relationships 

As  will  be  seen  below  the  ordinary  ambiguity  function  of  a  signal  such 
as  an  FM  or  CW  pulse  is  much  simpler  in  form  than  its  doppler  ambiguity 
function  while,  on  the  other  hand,  the  doppler  spreading  function  of  a 
target  is  much  simpler  than  the  ordinary  spreading  function.  Let  B  be 
the  bandwidth,  T  the  duration  and  fQ  the  carrier  of  the  signal. 

For  sufficiently  small  signal  time  bandwidth  product  BT  and  target 
time  stretch  a  sufficiently  near  1  ,  namely  for 

BT  | a  -  1|  <  0.2  (Eq.  6) 

the  frequency  stretch  a  may  be  approximated  by  a  frequency  shift 
(a-l)fQ  since  in  this  case  the  maximum  phase  error  caused  by  so  doing 

is  negligible. 

When  Eq.  6  holds,  it  follows  that 

$jq(x,v)  ~  j~  Sjq(t,  1  )  .  (Eq.  7) 

o  o 

Note  that  Eq.  7  is  exact  for  a  =  1  (i.e.  v  =  0),  since  the  approximation 
is  exact  for  a  =  1 . 

Clearly  Eq.  6  usually  holds  for  a  CW  signal  (with  BT  =  1)  in  typical  sonar 
and  radar  applications.  On  the  other  hand,  for  the  linear  FM  in  typical 
sonar  applications  as  opposed  to  typical  radar  applications  for  which 
| a- 1 1  is  usually  much  smaller  — (6)  often  does  not  hold. 
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From  the  foregoing  It  is  obvious  that  for  a  low  BT  signal,  such  as  a  CW 
pulse,  analysis  of  sonar  detection  and  measurement  extraction  may  be 
carried  out  using  the  ordinary  theory  with  the  exact  spreading  function 
replaced  by  an  approximate  spreading  function  found  from  the  Doppler 
spreading  function  of  the  target.  On  the  other  hand  for  a  high  BT  signal 
such  as  a  linear  FM  pulse  the  exact  spreading  function,  which  is  apt  to  be 
much  more  complicated,  must  be  used  if  the  ordinary  theory  is  applied.  One 
can  regain  a  simple  form  for  the  spreading  function,  at  the  cost  of  a 
moderate  increase  in  complexity  of  the  ambiguity  function,  by  using  the 
frequency  stretch  version  of  the  theory. 


There  is,  however,  one  fly  in  the  ointment  —  it  is  more  convenient  to 
represent  the  time  and  frequency  spreading  effects  of  the  medium  via  the 
ordinary  spreading  function.  Fortunately,  it  is  possible  to  approximate 
D  v  0  ■  2  f 

+  ft  f°r  v  < - as  an  ordinary  cross-ambiguity  function 

o  BT 

and  thus  show  that 


4>o(t,a) 


e27Tjv*  (T-T  '  ) 


~  <t> j(T-T',a-  7-)Stq(t',v')  dv'  dt * 

—  0 

f°  (Eq.  8) 


if  Sjq(t 1 ,v' )  is  negligible  for  v'  >  0.2  fQ/BT.  For  example, 

if  f  =  3,500  Hz,  B  =  220  Hz  and  T  =  0.5s;  this  limit  is  roughly 

7  Hz  (or  a  Doppler  shift  of  ±  3  knots),  which  should  be  sufficient 
to  include  most  frequency  spreading  in  the  ocean.  Use  of  this 
approximation  allows  us  to  use  the  Doppler-shift  version  for  large  BT 
signals  while  using  the  ordinary  spreading  function  to  model  the  medium. 
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2  SOME  APPLICATIONS  OF  THE  THEORY 

In  this  section  some  applications  of  the  theory  just  summarized  are  given 
both  to  illustrate  the  theory  and  because  of  their  own  inherent  interest. 
The  goal  of  these  applications  is  to  derive  the  relationship  between  the 
matched  filter  output  for  the  bistactic  echo  from  a  linear  FM  pulse  for  a 
moving,  turning  line  target  in  a  non-dispersive  medium.  To  apply  the 
fundamental  relationship  to  determine  the  matched  filter  output,  it  is 
necessary  first  to  use  the  definitions  to  obtain  the  ambiguity  function  of 
a  linear  FM  pulse  and  the  bistatic  spreading  function  of  a  moving,  turning 
line  target.  Along  the  way  it  is  convenient  to  include  the  propagation 
delays  in  the  target  spreading  function. 


2.1  The  ambiguity  functions  of  a  linear  FM  pulse 

As  an  example  of  an  ambiguity  function  consider  the  modulated  linear  FM 
pulse 


xFM^’^o’k’^ 


/T 


(t;T/2) 


2ttj  ( fQt+^kt2 ) 


(Eq.  9) 


where  f0  is  the  carrier  frequency,  T  is  the  pulse  length, 

B  =  | k j T  is  the  bandwidth  and 

J  L  (t;T)  =  j  1  if  N<T  .  (Eq.  10) 

!  0  if  j t j >T 

Its  ambiguity  functions  are 

lTi  2irj(f  +3*v)t 

yfm(t,v;  fo’k’T)  =  (1  '  •Lp)J"'(T.T)e  0  sinc[ir(T-|f| )(kx-v)] 

2ttJ  ( f  +’?v)x 

~  e  s i nc[xT  T(kx-v)]  for  x  «  T  . 


(Eq.  11a) 


2TTj(^)f0X 


U  / 

^FM  T,ot 


f0»k,T) 


'FM 


T  , 

[ax-(a-l)-|  ,  i  ;  fQ,k,T]  (Eq.  11b) 
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where  an  approximation  to  will  be  given  shortly.  The  ordinary 
ambiguity  yp^  given  in  Eq.  Ua  is  well  known  and  needs  no  comment  beyond 
the  fact  that  the  ordinary  ambiguity  function  of  a  CW  pulse  is  obtained  by 
setting  k  =  0  in  Eq.  11a.  The  doppler  ambiguity  function  yj^  is 
expressed  in  Eq.  lib  in  terms  of  y^  because  the  matched  filter  output 

for  a  linear  FM  signal  frequency  stretched  by  a  ,  from  Fig.  2  and 

D  ® 

Eq.  5b  (since  SjQ  in  this  case  is  6(x )  6(a-aQ)),  is 


yFM^ot’  a J 
0 


a  +1 

^j(~T—  )fr 


=  e 


'FM 


,  a  -1  fQ 
<T  Vr  ■  “o> 


(Eq.  12) 


Thus  if  (a  +1)  fjd.  is  taken  as  the  carrier  of  the  matched  filter  output 

M  o  o 

YpM  represents  the  modulation;  it  is  time  shifted  so  that  its  peak  in  x 
for  given  a  is  near  x  =  0. 

M 

When  YpjYj  is  found  by  substitution  of  Xp^  for  both  x  and  y  in  Eq.  lb 

and  use  of  the  result  in  Eq.  lz,  it  is  can  be  expressed  in  terms  of  Fresnel 
integrals  with  rather  complicated  limits;  however  suitable  approximation 
yields  simpler  limits  that  are  valid  for  most  sonar  purposes.  It  is 
convenient  to  give  the  results  in  terms  of  the  following  normalized 
variables: 

x  =  kT  /T7k0/a  *  kT  flJ?Q  , 
y  =  (oc- 1 )  /fj7SMa-l)  /y*  , 
v  =  a  SJJTJ*  /TTy*  , 

z  =  k  T  t  (Eq.  13) 


In  general  f  and  T  (and  hence  v)  are  determined  by  external 

circumstances  and  the  choice  of  k  is  available  for  optimizing  the 

properties  of  y^  ;  therefore  x  ,  y  and  z  are  referred  to  as  the 

normalized  bandwidth  (since  J kT |  =  B),  normalized  time  stretch  and 
normalized  time  shift  respectively. 

In  terms  of  these  variables 

-f  j* 

YpM(z,{/;  x,v)=kI~x[{y  -  T^)/x’1+^3-fi+^y7  xy 

*{c(b)-c(a)+j  sgn(x{/)[S(b)-S(a)}  +  e  (Eq.  14) 
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where 

C(t)  +  j  S(t)  = 


9 


_  -  oo 

a  =  /~\xy T  -  0+vy)  +  \  (y  -  TT^/)/x| 

b=  /  NT  jjj  +  d+vy)  -  \{y  -  )/x| 

*  =  £  -  4*  '  ^)/X  ’ 

e  <  |«/|  w  (y+^  +  0.27|x{/| ) 

and  sgn(  )  means  sign  of.  Not  that  for  v  =  0  these  relationships 
are  both  much  simpler  and  exact  (i.e.  e  =  0). 

It  will  be  seen  that  the  maximum  value  |a-l|  of  interest  is  0.02.  In 

terms  of  the  original  variables  e  <  |ot-l|(^+  0.27  BT  |a-l|):  there¬ 
fore  for  BT  <  750,  e  can  be  neglected  with  respect  to  the  maximum  of 

M  M 

YpM>  which  is  from  Eqs.  la,  lb,  11a  and  11b  YpM(0,0;  x,v)  =  1. 

Since  Yp^(aT,  1/a)  is  the  output  of  filter  matched  to  XpM  when  its  input 

is  Xpi^  frequency  stretched  by  a  ,  its  fourier  transform  with  request  to 

x  is  the  product  of  the  matched  filter  transfer  function  times  the 
fourier  transform  of  the  frequency  stretched  signal.  For  large  BT  -  which 
for  fixed  x  corresponds  to  small  v  -  these  both  have  the  same  form: 
quadratic  in  phase  and  rectangular  in  magnitude;  hence  the  fourier 

transform  of  Yp^(aT,  1/a)  with  respect  to  x  is  readily  obtained.  This 

result  may  then  be  time  and  frequency  shifted  in  accordance  with  Eq.  12  to 

yield  r^ff.a)  the  fourier  transform  of  Yp^(^,a)  with  respect  to  t. 

(If  the  fourier  transform  is  made  with  respect  to  the  normalized  variable 

z  instead  of  x  the  result  is  B  rtL  ;  it  is  convenient  to  ignore  this 

factor  and  consider  to  be  the  fourier  transform  of  Yp^  with  respect 

to  z  expressed  as  a  function  of  the  normalized  frequency  w  =  f/kT.) 

The  result  in  normalized  variables 


.(«W  *.*)*"<!(/*;  l)-<a»;  l-tfl)  e2"^  +  ''3/4x' 


This  expression  is  exact  for  v  -  0  and,  in  fact  its  inverse  fourier  with 
respect  to  u>  is  Eq.  14  with  .=  0  -  as  it  must  be. 
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To  determine  how  strong  the  dependence  of  and  are  on  v  , 

max  |YpMU0!/  ;  x,v)  -  YpM(z»t/  ;  x,o)|  was  computed  for  various  x  y 


and  v  .  The  results  are  shown  in  Fig.  4,  which  gives  this  maximum 
error  as  a  function  of  y  v/0.02  =  |a-l 1/0.02  for  various  x  and 
x/v  =  BT.  The  maximum  value  |a-l|  used  was  0.02,  since,  as  mentioned 
earlier,  this  will  be  seen  to  be  the  maximum  value  of  interest.  As 
might  be  expected  the  error  is  strongly  dependent  only  upon  BT. 


From  Fig.  4  it  is  clear  that  for  BT  >  50  the  error  in  using  v  =  0.0 

is  negligible  compared  to  the  peak  value  of  y^  ,  which  is 

Yp^O,  0  ;  x,v)  =  1.  One  should  not  expect,  however,  that  for  such 

values  of  BT,  Tp^  will  take  the  nice  form  given  in  Eq.  15.  In  fact 

because  of  Gibb's  like  phenomenon  the  actual  magnitude  and  phase  will 
oscillate  around  the  rectangular  magnitude  and  quadratic  phase  of 

Eq.  15.  Nevertheless  it  is  reasonable  to  say  that  rrM(w,«/  ;  x,o) 


approximates 
YpMU»</  ;  x,v). 


FMV 

M  M 

Ip|v|(w>!/>x5v)  whenever  yj!L(z,j/  ;  x,o)  approximates 


The  quadratic  phase  dependence  in  Eq.  15  implies  that  the  frequencies 
M 

in  YpM  differing  from  the  carrier  are  shifted  in  time  with  respect 
to  the  carrier;  i.e.  that  the  peak  of  y^  for  a  given  y  is 


broadened  in  y 


This  effect  is  shown  graphically  in  Figs.  5,  6  and  7, 
which  are  graphs  of  lYp^l2  versus  z  for  v  =  0  and  various  x  and 
y  .  In  fact  the  fourier  transform  of  1  ^ 


YppjU,^  ;  x,o)|2  with  respect 

y  is  readily  found  by  autocorrelation  of  FpM(«My  ;  x,o)  with 


to 

respect  to 


w  to  be 


?J-M{//x  ;1)J— L(w,  1  - 1  !//x  | )  (1  - 1  t//x  |  -  |w|  )sinc[2ux  y  w(  1  - 1  <//x|  -  |w| )]  . 

(Eq.  16) 


Two  limiting  cases  are  of  interest:  for  small  xy  the  last  factor 

is  essentially  unity  for  |w|<l-|x,t/|  and  |YpM|2  is  just  the  inverse 

fourier  transform  of  the  first  two  factors  with  respect  to  w  ;  namely 
( 1  - 1  f//  x  | ) 2  sinc2[u(l- \y/x\ )z]  -which  is  of  course  the  classical  result 
found  when  the  doppler  effect  may  be  approximated  by  a  frequency  shift. 
On  the  other  hand  for  large  xy  the  first  two  factors  in  Eq.  16  may  be 
replaced  by  unity  and  the  last  approximated  by  sin  [2ir|  xt/|w(l - 1 x/w| )]; 
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hence  ) Yp 


|t/x|(1-|(//x|)' 


[z  ;  %|f/x|(l-|t//x|)j  -  which  is  the 


result  obtained  by  Weston  [4], 


The  FM  doppler  ambiguity  function  just  discussed  in  detail  is  for  a  pulse 
whose  center  corresponds  to  t  =  0.  In  section  2.4  the  FM  doppler 

ambiguity  function  Yp^(ctT,  1/a)  for  a  pulse  whose  center  has  an  abitrary 
time  ty  is  required.  Furthermore,  in  that  section  a  simple  relationship 
between  the  doppler  ambiguity  function  for  a  close  to  aQ  and  that  for 
aQ  is  required.  If  the  definition  of  the  doppler  ambiguity  function 
(Eq.  lb)  is  applied  for  an  FM  pulse  centered  on  ty  and  the  frequency 
stretch  a  replaced  by  a  frequency  stretch  aQ  plus  a  frequency  shift 
(a-aQ)f0  ,  the  result  after  suitable  algebric  manipulation  and  approx¬ 
imation  is 


D 1 

YFM(ax,  1/a) 


a-a„  f  2 

e  o 


WT  +  a2  k  fo  +  a  V  ’  VaQ] 


which  is  valid  if 

i  i  .0.2  o.2B  nr 

|a-aQ I  <  mi n  gj-  ,  0.05  ^  ^  /  gy 


(Eq.  17) 


(Eq.  18) 


Factorinc 


;ion  delays  into  the  target  spreading  function 


An  important  application  of  the  fundamental  relationships  illustrated  in 
Fig.  3  is  to  make  a  convenient  modification  to  the  medium  and  target 
spreading  functions.  In  their  most  straightforward  versions  the  propaga¬ 
tion  delays  to  and  from  the  target  are  included  in  the  medium  spreading 
functions.  It  is  more  convenient  to  transfer  the  direct  path  delays  into 
the  target  spreading  function  so  that  with  a  non-dispersive  medium  only  the 
target  spreading  function  need  be  used  to  relate  the  matched-filter-bank 
outputs  to  the  ambiguity  function  of  the  transmitted  signal.  In  this 
version  the  medium  spreading  functions  are  only  used  to  represent  dispersive 
effects  of  the  medium. 


The  procedure  to  be  carried  out  is  illustrated  in  Fig.  8  for  the  frequency 
shift  version  where  Sq^  , ,  S^,  and  are  the  spreading  functions  for 


FIG.  8  RELATIONSHIPS  BETWEEN  TWO  VERSIONS  OF  A  SIMPLIFIED  ACTIVE  SONAR  MODEL 
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the  straightforward  version  (delays  in  the  medium  spreading  functions); 

S01*  S12  and  S23  are  the  sPreadin9  Unctions  for  the  more  convenient 

version  (delays  in  the  target  spreading  function)  and  SjD(t,v;  Td)  and 

Sj^(t,v;  id)  are  the  spreading  functions  of  a  time  delay  and  a  time 

advance  by  xd.  Obviously  Sy^  followed  by  SyD  or  vice  versa  cancel 

each  other.  The  results  from  use  of  Fig.  3  and  Eq.  5a  to  combine  spreading 
functions  are  as  follows: 

SqI^^)  =  1  (T+T(jj»v)  > 

S]2(t,v)  =  e"2lTjTdT  $12(x-TdT-TdE,v) 

S23(t,v)  =  e27r^TdE  S2>3(T+xdE,v)  ,  (Eq.  19a) 


where  xdy  and  xdE  are  the  transmission  and  echo  direct  path  propagation 
delays.  From  Fig.  8  it  is  obvious  that  S12(x,v)  is  the  spreading 
function  of  the  overall  system  for  a  non-dispersive  medium. 

For  the  time  stretch  version  the  same  procedure  may  be  used  to  yield 


Sqi (t ,a)  soi 1 (T+TdT»a)  * 
sj2(x,a)  =  1 2 .  (t — jjp- -  TDE,a) 

S23(T,a)  =  S2,3(x+  — p,a)  , 


(Eq.  19b) 


as  before,  S^2(x,a)  is  the  Doppler  spreading  function  of  the  overall 
system  for  a  non-dispersive  medium. 


2.3  The  bistatic  spreading  function  of  a  moving,  turning  line  target 

The  target  model  and  geometry  to  be  investigated  is  shown  in  Fig.  9,  where 
the  variables  used  are  defined  in  Table  1.  The  values  of  the  variables  of 
interest  are  those  that  obtain  at  a  given  time  that  may  conveniently  be 
taken  to  be  t  =  0.  It  is  assumed  that  the  medium  is  noise  and  reverbera¬ 
tion  free,  has  no  propagation  loss  and  does  not  spread  the  sianal  in  time 
or  frequency.  Taking  into  account  noise,  reverberation  and  propagation 
loss  presents  no  particular  problem  and  may  be  done  essentially  independently 
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FIG.  9  BISTATIC  GEOMETRY 


TABLE  1 


LIST  OF  VARIABLES 


a 

Target  aspect 

Y 

Bistatic  angle 

*T 

Bearing  of  target  with 

respect  to  transmitter 

Bearing  of  target  with 

respect  to  receiver 

V 

Target  speed 

VT 

Transmitter  speed 

VR 

Receiver  speed 

rT 

Range  from  transmitter 

to  centre  of  target 

rR 

Range  from  receiver  to 

centre  of  target 

a 

Distance  along  target 

PT(o) 

Range  from  transmitter 

to  point  at  o  along  target 

PR(o) 

Range  from  receiver  to 

point  a  along  target 
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of  the  development  presented  herein;  however  including  time  and 
frequency  spreading  requires  a  modification  of  the  present  development 
using  the  techniques  described  above.  One  point  of  possible  confusion 
is  the  use  of  a  for  both  the  aspect  and  the  doppler  shift;  fortu¬ 
nately  the  context  makes  clear  which  use  is  implied. 


For  short  periods  of  time  each  point  along  the  target  may  be  assumed 
to  have  constant  range  rates  with  respect  to  the  transmitter  and 
receiver  that  are  readily  determined  from  the  geometry.  The  corres¬ 
ponding  time  stretches  and  hence  the  doppler  spreading  functions  for 
each  point  may  then  be  determined  and  summed  to  yield  the  target 
spreading  function.  The  result  under  the  assumation  that  the  signal 
arrives  at  the  target  near  t  =  0  is  for  a  target  whose  points  have 
reflectivity  f (a) 


SD(x,a) 


f-cx+l-—rr  \ 


a£l 

a 


(Eq.  20) 


where  6(  )  is  the  Dirac  delta  function,  c  is  the  speed  of  sound  and 
1  *  rT  +  rR  ’ 

i  -  rT  +  rR  =  -( v<5  +  vT  cos  ipT  +  vR  ipR)  , 

6  -  [cos  a  +  cos(tx-y)]  , 
a  -  i  cos  a/cS  , 

6  =  [a  sin  a  *■  (a-y)  sm(a-y)]/6  .  (Eq.  21) 


For  a  monostatic  situation  the  same  results  apply  with  y  :  0  , 

rT  =  rR  =  r  ,  t^y  =  ^R  =  ^  anci  v’i  =  vr  =  VTR  »  hence 

l  2r  , 

i  -  2r  =  -2(v  cos  a  +  v^R  cos  'p)  , 

6  =  2  cos  a  , 
a  =  if  2c  , 

6  =  a  tan  a  ,  (Eq,  22) 
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2.4  The  matched  filter  output  for  the  bi static  echo  from  a 
linear  FM  pulse  for  a  moving,  turning  line  target  in  a 
non-dispersive  medium 

To  find  the  match  filter  output  d>^(x ,1 )  it  is  necessary  to  substitute 

the  approximate  FM  ambiguity  and  target  spreading  functions  giver,  in 
Eqs.  17  and  20  into  the  fundamental  relationship  defined  by  Fig.  2  and 
Eq.  5b.  Recall  that  ty  is  the  center  of  the  FM  pulse.  If  ty  =  -ry/c 

the  signal  will  arrive  at  the  centre  of  the  target  at  t  =  0  ;  therefore 
let  Aty  =  ty  +  ry/c  so  that  Aty  =  0  implies  ty  =  -ry/c  .  The 

results  to  be  given  only  apply  for  small  Aty  say  on  the  order  of  the 

length  of  a  single  ping  (which  may  include  several  pulses)*  since  only 
then  will  (Eq.  20)  be  valid.  The  result  of  the  substitutions  with 
a  =  1  -  l/c  is 


n  f-(x-Tj)  M 

<}>y(T,l)  %  e  iAt-t^) 


(Eq.  23) 


where 


U  +  A  ty  )  /  C  , 


Td  »  U  +  ^r(Aty  +  ^)]/c  , 

4%)  =  r  [At)]  ■  rl(f,  i-i/c)HM(f)  , 


rfM(f,a)  =r[  Y”M(T,a)]  , 


HM(f)  =  F ( S { [  1  -  AtT  -  a](f+f ')  -  —  (f  +  f_))/c)  , 

(c-*)*  T  0  (c -l)2  k  2c  0 


F(f)  =  T  [f(-T)]  , 


(Eq.  24) 


and  T  []  indicates  Fourier  transform  of  . 


In  general  t ^  may  not  be  known.  Tkii  it,  not  impofitant  60  long  ai 
the  iignal  and  matched  ^ilten  one  toned  to  the  iame  *e(,eAence. 
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Several  1  points  may  be  observed  about  this  hodgepodge: 

.  The  carrier  of  the  matched  filter  output  is  both  time 
and  frequency  shifted  as  a  function  of  l 

.  The  modulation  of  the  matched  filter  output  is  shifted 
in  time  by  a  different  amount  as  a  function  of  i 

.  Transformed  into  the  frequency  domain  the  modulation 
of  the  matched  filter  output  is  the  frequency  shifted 
and  stretched  target  transfer  function  F(f)  observed 

through  the  window  Fp^f),  which  has  been  discribed 

above. 

.  Since  [c/(c+£)]2  %  1,  the  effect  of  3  on  the  transfer 
F(f)  is  a  stretch  in  frequency  proportional  to 
l/[i-3(At+fo/k)]  %  l+3(At+fQ/k)  and  a  shift  in  frequency 

of  3(Atp  fg  +  £  f2/7ck)*gAtp  fQ  since  |£|/c  <  0.02. 

At  this  point  a  word  about  the  validity  of  the  approximation  used  above  is 
in  order.  In  the  first  place  Eq.  18  must  be  obeyed  if  Eq.  17  is  to  be 
valid  and  hence  Eqs.  23  and  24  hold.  If  the  length  of  the  target  is  L 
then  f(a)  must  be  zero  for  o  >|L|;  therefore  from  Eq.  20  the  maximum 

value  of  |a-aQ|  is  (recall  aQ  =  l-£/c)  |8  lj-6(l-a)|  ,  whose  maximum 

in  turn  is  L[o|.  Thus  for  fQ  =  3500  Hz,  6  =  220  Hz,  T  =  0.5s  and 

L  =  100  m,  Eq.  17  implies  |a|  <  3.82  /s  corresponding  to  a  full  turn  in 

1.57  minutes.  If  Eq.  6  is  satisfied  then  the  ordinary  ambiguity  function 

may  be  used  to  determine  IpM  .  The  value  of  a  of  interest  in  Eq.  6 

is  a  =  aQ  ;  hence  a-1  =  -£/ c  ,  whose  maximum  is  (2v  +  Vy  +  vR)/c. 

Thus  for  the  signal  parameters  just  given  |2v  +  Vy  +  vR|  <  2.73  kn  to 

Insure  the  ordinary  ambiguity  function  may  be  always  used  -  a  rather 
severe  constraint.  On  the  other  hand  constraining  |2v  +  Vy  +  v|  to  be 

less  than  60  kn,  which  is  equivalent  to  ja-1 |  <  0.02  ,  does  not  seem 

very  onerous.  This  constraint  has  been  used  previously  in  discussing 

approximations  for  Yp^  • 
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CONCLUSIONS 


Two  version  of  linear  time-varying  system  theory  are  of  use  in  analysis  of 
active  sonar  systems:  the  ordinary  (frequency-shift)  version,  which 
applies  for  CW  signals  and  the  doppler  (frequency-stretch)  version,  which 
applies  for  FM  signals.  When  the  doppler  version  is  used  to  determine  the 
matched  filter  output  for  the  bi static  echo  from  a  moving,  turning  line 
target  in  a  non-dispersive  medium  the  result  is  that  the  modulation  of  the 
matched  filter  output  is  the  time-shifted  convolution  of  a  frequency 
shifted  and  stretched  impulse  response  representing  the  target  and  a  window 
function  representing  the  signal. 
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APPENDIX  A 

DERIVATION  OF  THE  RELATIONSHIP 
BETWEEN  THE  IMPULSE  RESPONSE  AND  THE  SPREADING  FUNCTIONS 


Substitution  of  yj(t)  =  into  Eq.  2a  yields 


h(t.t-T')  =  //e2"J''(t'T,6(t-T-T')S(v,T) 


-00-00 

since  h(t,t-x')  is  the  output  due  to  an  impulsive  input  at 
t-(t-x')  =  x'.  Replacing  t-x1  by  x(=>  x'  =  t-x)  yields  Eq.  3a. 
Similarly  substitution  of  yx ( t)  =  6( t-x 1 )  into  Eq.  2b  yields 

00  00  _ 

h,n(t,t-x')  =  f  f  6[a(t-x)  -x * 3  Stn(T»a)  dadx 


c  l-L  a  6{oit  'T'T,)  SIO(a’a)  dtda 

_  7  1  cD  /ott-x1  \  j 
-  o  J  SI(T~a  ’  «>  dr‘ 


Replacing  t-x‘  by  x  yields  Eq.  3b. 
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APPENDIX  B 

DERIVATION  OF  THE  FREQUENCY-SHIFT  VERSION 
OF  THE  FUNDAMENTAL  RELATIONSHIPS 


The  derivation  takes  place  in  two  steps:  first  some  properties  of  the 

modified  double  convolution  (x) 1  are  derived  and  then  these  properties 

x.v 

are  utilized  to  prove  the  relationships. 

B.1  PROPERTIES  OF  THE  MODIFIED  DOUBLE  CONVOLUTION 
From  the  definition  Eq.  5a 

a  (x) '  b(x,v)  =  /  /  e^v  *  ^b(x-x '  ,v-v' )a(x 1 ,v‘ )dv'dx 1 

x,v  _0° 

=  /  /eilT^v  v  ^  a(x-x' ,v-v')b(x' ,v')dv'dx' , 

-00-00 

which  clearly  does  not  in  general  equal  b  (x) '  a(x,v).  Furthermore, 

t,v 

from  Eq.  5a  with  x,v,x'  and  v'  replaced  by  x-x ' ,v-v' ,x"-x '  and 
v"-v' : 

b  (x)  a(x-x '  ,v-  v')  =  /  /  e^j(v"-v' )  (T-T"  ^a(x-x",v-v")b(x"-x  '  ,v"-v'  ,'dv"dx" 

X  , V  -oo-oo 

Substitution  of  this  into  Eq.  5a  yields: 

00  oo  00  00 

C  ®'  (b  X'  a)!x,v)  =  S  1  e2Tjv'(T-T>)  ,  ,  e2»j(w"-v‘)(T-T*J 

—on  *“00  —  00  “Of' 


x,v  x,v 


•  a(x-x" ,v-v")b(x"-x,v"-v' )dv"dx"  c(x 1  ,v’ )dv'dx  ' 

=  7  7  eZ,jv"(T-Tl'>a<T-T\v-v") 

-00  -00 

'  -h-L  e^'V  T  )b(x"-x!,v"-v,)c(x,,vl)dv,dx,dv,,dx" 

=  (c  (x) 1  b)  (x)  1  a(x,v)  , 
x,v  x,v 


on  use  of  Eq.  5a  twice  more. 


**>»«  wti 
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Thus  the  operator  ®  has  two  properties  of  importance:  it  is  associative, 
but  not  cummutativeI’vThe  true  double  convolution  operator  ®  ,  which 

T,V 

differs  from  ®'  in  that  it  is  defined  by  Eq.  5a  with  the  exponential 

T,V 

term  of  the  integrand  removed,  is  both  associative  and  commutative. 


B.2  THE  FUNDAMENTAL  RELATIONSHIPS 

Note  the  similarity  in  form  between  Eqs.  5a  and  2a,  the  definition  of  the 
spreading  function.  The  definition  of  the  cross-ambiguity  function  (Eq.  la) 
is  seen  to  take  the  same  form  if  it  is  rewritten 


*l(x,v)  =  _/_/  e^v  T  ^ [6 (v-v ' )  x*(T-T,)]yi(T)dv,dT*. 


As  a  direct  result  both  Eqs.  la  and  2a  may  be  written  in  terms  of  ®  : 

t,v 

<^(t,v)  =  y-  ®  '  x(t.v)  , 

T,V 


y0(x)  =  SI0  ®'  yj(x,v)  ,  Eq  B1 ) 

x,v 


where 

x(x,v)  =  x*(x)6(v) . 


Use  of  these  two  results  and  the  associative  law  yields: 
<{>0(t.v)  =  y0  x  It.v) 

T,V 


=  (SI0  ®  yj)  ®  x(x.v) 
x,v  x,v 


=  SI0  &  ®  x)(t,v) 

x,v  x,v 


S jq  ®l  41  j  ( t  , v) , 

X,v 
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which  is  the  relationship  of  Fig.  2.  Use  of  Eq.  B1  twice  and  the 
associative  law  yields: 


ycM  -  ^bc  ®  ^ab  ®  yA)(T»v) 

T  ,V  T ,V 

=  ^BC®  ^AB^  ® 

T,V  T,V 


But  from  Eq.  Bl, 


yC(r)  =  SAC 


yA(T.v). 

T,V 


Since  these  results  hold  no  matter  what  y.(x)  is,  it  follows  that 
the  relationship  of  Fig.  3  must  hold. 
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APPENDIX  C 

DERIVATION  OF  THE  FREQUENCY  STRETCH  VERSION 
OF  THE  FUNDAMENTAL  RELATIONSHIPS 


As  with  the  frequency-shift  version,  the  derivation  takes  place  in  two 
steps:  first  some  properties  of  0  are  derived  and  then  these 

x,a 

properties  are  utilized  to  prove  the  relationships. 


C.l  PROPERTIES  OF  THE  MODIFIED  DOUBLE  CONVOLUTION 

From  the  definition  Eq.  5b  and  an  interchange  of  order  of  integration: 


aD  ®'  bD(x,a)  -  7/°bD[a'(-i-x'),  £,]  a^x'.a'Jdx'  da' 

n  -oo  ot 

T,a 


o  -°° 


-  7  7'-.  aD(T bD(T',a')dt'  da'  , 

0  -OO a  a  u 

which  clearly  does  not  in  general  equal  bD®  aD(x,a).  Furthermore, 

T,a  a 

from  Eq.  5b  with  T,a,x'  and  a!  replaced  by  a1 (x-x ' ) ,  ,a* (x"-x ' ) 

and  rr  : 


OO  OO  -  I! 

bD®  aD[a'  (t-x * ) ,  2fr]  =  /  /  aV(x-x“),  J  bULct '  (t"-t ' ) ,  jMda"  dx" 

U  -oo  o  U  U 


T,a 


PKSCEDING  PAS®  BUNK*M<yt  n 
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Substitution  of  this  into  Eq.  5b  yields: 


n  n  oo  oo  co  do  n  H 

cD  ®  (bD®  aD)(x,a)  =  ffff  aV(T-T"),  V)bV (t"-x'  ) ,  V|da"  dx" 

-OO  0-00  O  U-  Ub 

T,a  i,a 

.  c^(Tl,a')dal  dx’ 

co  co  p.  00  00  n  11 

=  f  f  aV(T-i"),  ~rr]  f  f  bD[a  *  (x" -t  * ) ,  V| 

-oo  o  u  -oo  o  a 

.  c^(t' ,a' )da'dx'da"dT" 

=  (cD  ®  bD)  (x)’  aD(x,a)  , 

T,a  x,a 

on  use  of  (5b)  twice  more.  Again  x  1  is  associative  but  not  commutative. 

x  ,a 


C.2  THE  FUNDAMENTAL  RELATIONSHIP 


Note  the  similarity  in  form  between  (5b)  and  (2b),  the  definition  of  the 
doppler  spreading  function.  The  definition  of  the  doppler  cross-ambiguity 
function  (lb)  is  seen  to  take  the  same  form  if  it  is  rewritten 


rs  00  00  | 

<J>^(t,oi)  =  f  f  5{~~)  x*[t-t')3  y.(x' )  da*  dx' 


(The  factor  a  In  front  of  the  integral  in  the  original  definition  (lb)  was 
inserted  to  yield  this  convenient  form.)  As  a  direct  result  both  (lb)  and  (2b) 
may  be  written  in  tem.s  of  ® 1 : 

i,a 

><*)  =  yj  ®  '  XD(T,a), 
x,a 

y|5u)  =  sio  ®  '  yr(T,a),  (Eq.  Cl) 

x  »a 

where 

XD(x,a)  =  x*(x)6(^— ) . 
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Use  of  these  two  results  and  the  associative  law  yields: 

<J>[|(T,a)  =  yQ  x  '  x  D(T,a) 

T,a 


■  (Sj5  ®'  y.)  ®'  xD(t,a) 
x,a  i,a 

=  S°  ®'(yj  ®  1  xD)(T,a) 
x,a  x,a 


=  S  D  ®  *  <j>i  (t >a)  , 

10  T,a 

which  is  the  relationship  of  Fig.,  Z.  Use  of  (B-l)  twice  and  the  associative 
law  yields: 


yc(x)  =  S  D  ®  '(SA°  ®  *  yA)  (r.a) 
BC  T,a  x,a 


(SBC  ®  '  SAB)  ®  1  yA(T,a) 

T,a  x,a 


But  also  from  (C-l), 


yc(T)  =  sAr  ®  '  yA(x.a). 
x,a 


Since  these  results  hold  no  matter  what  yA(x)  is,  it  follows  that  the 
relationship  of  Fig.  3  must  hold.  H 
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APPENDIX  D 


DERIVATION  OF  THE  PPROXIMATE  RELATIONSHIPS 


This  appendix  is  concerned  with  showing  two  things:  For  low  3T  product 
signals,  such  as  CW  pulse,  the  spreading  function  of  a  moving  target  may 
be  replaced  by  its  much  simpler  doppler  spreading  function  and  the 
standard  version  of  the  theory  applied.  For  high  BT  product  signals, 
such  as  linear  FM  pulses,  the  doppler  spreading  functions  of  the  medium 
may  be  replaced  by  its  ordinary  spreading  functions  and  the  doppler  version 
of  the  theory  applied. 


D.l  REPLACEMENT  OF  THE  ORDINARY  SPREADING  FUNCTION  BY  THE  DOPPLER 
SPREADING  FUNCTION 


A  fundamental  approximation  used  in  communications  theory  is  that  for  a 
small  BT  signal  modulated  on  a  carrier  x(t)  the  doppler  shift  affects 
only  the  carrier: 

2irj(a-l  )f  t 

x (at )  *  e  0  x(t)  .  (Eq.  D1 ) 


I 


This  approximation  is  discussed  in  detail  by  Rihaczek  {[i]  pp  56-65),  where 
he  derives  the  general  condition  that  for  it  to  be  a  good  approximation  (6) 
must  be  satisfied. 

When  the  approximation  is  valid  its  substitution  into  (lb)  yields 

00 

n  f  -2irj (a-1 )  f  (t-x) 

<^(x,a)  x  a  I  e  x*(t-x)  y^t)  dt  , 

-00 

which  on  comparison  with  (la)  becomes 
<J>!?(T»a)  a4>i [t , (a-1 )  fQJ. 


Substitution  of  this  result  into  the  fundamental  relationship  given  by  Fig.  z 
and  (5b)  yields: 


¥x’  (a'1J  fo]  *  a!]  «•  4)I[a,(T"T,)’  ^'"1)fo]  dx’ 


oo  oo  2tt j  (a‘  -1  )f  (t-t  ' ) 


/  /  e 

-00  0 


4>j['f-x 1 ,  (a-a1  ^ f 0 3S jq ( x ’  .a1  )da'dx' 


BI/NK.NOT  FI- 
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Since  from  (la)  and  (D-l) 

<j>T(a'x,v)  =  f  e"2lTjv(t‘a,T)  x*(t-a'x)  y.(t)  dt 
1  -00  * 

=  a-  /e"2lTjva'(t‘T)  x*[a'lt-x)]  y,(a't)  dt 

-CO  I 

2nj(a'-l)f0T  .  2.ja'v(t-t) 


a’e  '  /e 

-00 

2irj(a'-l)f  ■ t 


x*(t-x)  yj(t)  dt 


a'  e 


4>  j  (x  ,a'v) 


Comparison  of  the  above  equation  for  4>q[x , (ct-1 ) f Q]  with  the  fundamental 
relationship  given  by  Fig.  2  and  (5a)  yields  (7). 

D.2  REPLACEMENT  OF  THE  DOPPLER  SPREADING  FUNCTION  BY  THE  ORDINARY 
SPREADING  FUNCTION 

An  obvious  extension  of  D-l  is 

2frj  (a'-a)f  t  n? 

x(a't)  «  e  0  x(at),  |a'-a|<  .  (Eq.  D2) 

To  derive  (8)  from  this  result  note  its  substitution  into  (lb) 

— 4r  *i(t.a+  r>*"  e‘!”jv(t-T)  x*[a(t-t) J  y,(t)  dt 

a  +  |-  o  _o° 

o 

for  v  <  0.2  fQ/BT.  The  right  side  of  this  equation  is  the  ordinary  cross¬ 
ambiguity  function  between  x(at)  and  yj(t).  Furthermore,  from  (lb), 

<j>q(x,cx)  is  the  ordinary  cross-ambiguity  function  between  x(at)  and 

yQ(t)  evaluated  at  v  =  0;  therefore  substitution  of  <j>q(x,cx)  for 

<Pq(t ,0)  and  •  <J>j(x-x '  ,a-^-)  for  «i>j (x-x '  ,-v 1 )  in  the  fundamental 


relationship  given  in  Fig.  2  and  (5a)  with  v  =  0  is  permitted  and  yields  (8). 
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APPENDIX  E 

DERIVATION  OF  APPROXIMATE  FORMS  OF  THE  LINEAR  FM  DOPPLER  AMBIGUITY  FUNCTION 


Three  approximations  are  to  be  derived  in  this  appendix:  The  first  two  are 
approximations  to  the  window  function  y^  from  which  the  FM  doppler 
ambiguity  function  is  readily  obtainable  via  (lib).  The  last  approximation 
is  of  Yp|^(cxT,  1/a),  the  doppler  ambiguity  function  of  an  FM  pulse 
centered  at  ty  rather  than  0,  in  terms  of  l/aQ)  for  a  near  a  . 


E.l  THE  WINDOW  FUNCTION 

From  the  definition  of  the  doppler  ambiguity  (lb) 

°° 

YpM(ai,  1/a)  =  1  /  x*  -x  )  x  (t)  dt 
m  FM  a  FM 


CO 

=  /  X*  (t-T+x')  X  foe  ( t  ♦*!  ’  >3  dt 
FM  FM 


after  change  of  dummy  variable.  Substitution  of  this  result  and  the  definition 


of  the  FM  pulse 

Xp^  into 

(12)  (with  aQ 

replaced  by  a)  yields 

M  ,  ,  a+1 

yfn  (t+  — 

v>- 

T  u(tMAX_tMlN^ 

^MAX  2ttj> 

e  ‘dt 

lMIN 

where 

u(t)  -  I  0 

t  <  0 

)  1 

t  >  0 

<j>0  =  ~f Q ( t~x+x 1 )  -  j  k(t-x+x')2+  afQ(t+x')  +  j  ka2(t+x’) 


a+1 

2 


f  i 

o 


=  (a-1)  fQ  (t-^+x 1 )+??  kai(t+x’ )?  -  j  k(t-x+x ' )2 
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and,  since  Xp^  (t-T+t')  is  zero  outside  the  interval 

and  Xp^[a(t+T' )]  is  zero  outside  the  interval  -t'+ 


tMAX  =  min[T'T' 


* 


k>  ‘i T  -  ^ 


*MIN 


max[t- 


T  -77 


T 

2’ 


liadi 

2a 


-]  • 


If  x1  is  chosen  to  eliminate  the  terms  of  <b  that  are  of  first  order 

in  t  then  can  be  easily  written  in  terms  of  Fresnel's  integrals. 

Rather  than  proceed  directly  on  this  track  it  is  appropriate  to  make  some 
approximations. 


Consider  first  the  following  approximate  limits: 


t' 


MAX 


) 


t' 

MIN 


The  maximum  error  between  the  actual  and  approximate  limits  is  T|a-1  |/a  ; 
furthermore  the  integrand  of  (E-l)  has  unity  magnitude;  therefore,  since 
u ( t ' MAX  '  t *min^  =  j-x(t’T)’  (E-U  may  be  rewritten 


M  /T+  a+l^SL  a)  s  Ltl 
yFMu  a  k  ,a  “ 


(t,T) 


,t’ 


t' 


MAX  +2irj(j) 


MIN 


dt  +  e 


where 

I ci  |  <  j  •  1  =  |a-l  |/a  =  v|y|/(l+v  y) 

from  (15). 


The  next  approximation  is  to  neglect  (a-1 ) 2  terms  ir.  <J>  ;  i.e.  to 
approximate  a2  =  (a-1 +1 ) 2  =  (a-1)2  +  2(a-l)  +  1  by  2(a-l)  +  1  =  2  a-1. 

If  <j/  is  <pQ  with  this  approximation,  clearly 

<P0  =  +  k  k(a-l)2  (t+  t')2  . 

Furthermore,  since  |eJX-l|<x  for  real  x  ,  (E-l)  may  rewritten 


M  /  a+]_  _o  \  _  1 

WT  a  k  ,c^  ~  T 


(t;T) 


t1 


dt  +  Cl  +  £  2 


MIN 
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where 


M  $  rJ“L(T;T) 


tMAX 


t1 
MIN 


tt | k {  (a-l)2(t+x')2 


The  value  of  this  last  integral  is  independent  of  t '  ;  it  is  convenient 
to  select  t'  -  x  for  its  evaluation  : 


|e  i  *  f  jk|  (a-1)2  t: 


T-|x|+x 


T- I T | -T 


=  ^(t;T)  f  |  k  j  (a-1 ) 2  (£)3+  (~-\i\)3 

'  n  IU  T2  (a-1 )2<0.27  BT(a-l)2  =  0.27  v  x  y2 


12 


from  (15).  Therefore  for  e  =  ex  +  t2 


|c|  <  |c,|  +  \C2 1  <  +  0.27  BT  |a-l|)  -  vlyK^t  O.Z7|x  y|)  , 

as  is  given  in  (16). 


Now  x'  is  selected  so  that  the  multiplier  of  t  in  <J>  is  zero;  i.e. 
so  that 

W' 

-KX-  I  =  (a-1)  f  +  ( 2a-l )  k  i'  -  k(x'-T)  =  0  .  (Eq.  E2) 

t=0 

If  this  is  multiplied  by  (t  -  j  +  x')  and  subtracted  from  <j>'  ,  the 
result  is  L  0 

^  =  («•  1)  k  t2  +  (t ' ) 2 -t ' (2!t ' -t ) ]  -  kL(x'-x)2-(x'-x)(2x’-x)  J 

=  (a-l)k  t2  -  (a-l)k  x'(x'-x)  . 

But  from  Eq.  E.2 

2(a-l)k  (x1-  j)  =  akx  +  (a-1)  fQ  ; 
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therefore 


t'0  =  (a-1)  k  t2- 


[akx+  (a-1 )f  j2  (a-1) 


4(a-l )k 


,  f  _ _ 

and  (E.l)  becomes  (on  replacing  t  +  — -  ^  by  x  and  2/  |  (a-1  )k|  t  by  t) 

it  . ,  b  , 


( I  a )  =  --  -r_L  (  t  -  2Ll3 — 2.  t  )  _§ — 
Yfm  (T,a;  T  (t  a  k  '  /rj 


J4> 


24(^1  )k| 


2  rJ  si n[ (a-1 )k]t 


dt  +  c 


where 


b  =  /  (a-l)k 


akx  T  ,  |T 
(a-1 )k  ^ 


-*£&]  ■ 


This  result  is  equivalent  to  Eq.  16  on  substitution  of  Eq.  15, 

When  BT  is  large  an  even  simpler  approximation  may  be  found.  As  noted 

in  the  main  text  the  Fourier  transform  of  Yp^ax.l/a)  is  the  product 

of  the  matched  filter  transfer  function  and  the  Fourier  transform  of  a 
frequency  stretched  linear  FM  pulse.  But  for  large  BT  (see  Rihaczek 
[1],  Chapter  7)  the  frequency  stretched  linear  FM  has  Fourier  transform 

S7T  "  <f  -  “fo  •  “B«>  « 


Furthermore,  since  the  matched  filter  is  matched  to  the  unstretched 
signal  it  has  transfer  function 

»j[j“(f-f0)  -  sin(k)/4] 

—  •r"Mf-f_  ;B/2)e 

/T 
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Taking  their  product  yields 


rFM<f’  W  =  ^ 


7TJ 


g-1  f  .a+1 
2  o  5  4 


B 


“Id  f2  _  f  f 

9  a  o 


m  n 

To  obtain  ,  it  is  apparent  from  Eq.  12  that  PpM(f,  1/a)  must 

be  frequency  shifted  by  -  fQ  and  then  time  shifted  by  -  —  . 

Since  in  the  frequency  domain  a  time  shift  of  Ax  is  a  multiplication 

by  e^nJ  ,  the  result  is 


1  o-l  f  . 

aB  2  o  * 

a+l  R 

~r B 

f  - 

la^iL  (.«h  Q . 

_^(a2-1) 

e 

V 

*{(0.-1)^ 

1  4 

,2+y^ 

w  k/i  .  kvy. 

(l+vy)B  x 

’  2 

Cri  ' 

v  2x  •  nr 

Trj(2+vy)xy  w2 

a 

vw.y2 

(2+vy)x 

*1 

i  > 

X 

on  use  of  (15).  Setting  v  =  0  yields  (17)  except  for  the  factor  1/B, 
whose  presence  was  explained  earlier. 

M 

The  Fourier  transform  of  |ypM(z,y;x,o) |  with  respect  to  z  is  from 

Fourier  transform  theory  just  Tp*(-w,y;x,o)  x  rpM(w,y;x,o)  where  x 

w  w 

signifies  convolution  with  respect  to  w  ;  i.e.  it  is  the  correlation 
of  TpM  with  itself.  Use  of  (17)  yields 

M  M  °°  M*  M 

r FM ^ ~w »-y » x » 0 ^  Pp|^  (w  y;x,o)  rp|j(w+w'  >y\x,o)  dw1  , 
w 

-00  *  X  * 
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If  0  <  W  <  1  -  |^|  , 

(2w';l-|f|)  [2(w+w');l-  |J|]  =  u(2w'+1-|J|)-u[2(w,+w)-l-|^|] 

and  the  integral  becomes 


(l-|J|-w)/2 

(  27Tjxy[(w+w')2-(w')2J  .  2iTjxyw2  e4TT^xyww 
J  e  aw  -  e  47rjXyw 

-O-ljD/2 


(l-|£|-2w)/2 


-H-ljfU/2 


2tt  jxy  [  ( 1  - |£| )  w-w2  ]  -2ttj  xy[  ( 1  - 1  )  w-w 2  ] 

i _ x  e _ x 

4irjxyw 


=  (l-ljl-w)  sin  2irxy[(l-|y-| )w-w2J  , 

which  agrees  with  (18).  The  same  procddure  may  be  carried  out  for 
-l+|y|  <  w  <  1  to  obtain  the  complete  result. 

A 


E . 2  THE  AMBIGUITY  FUNCTION  FOR  TIME  SHIFTED  LINEAR  FM  PULSES 
From  (D-l)  for  |a-a  |  <  0.2/BT 

xFM(at)  *  e2lJ  XpM(aQt)  , 


since  the  carrier  of  Xp^  is  frt  .  On  the  other  hand,  by  the  definition 


XFM’  ^ 


’WV  *  e 

0 


[a -a  1 

— r~  f2  +  (a-a  )f  t 
a  k  o  '  o'  o  J 


xFM(%t) 
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a -a 

provided  that  (a  t+ — r^-  f  ;T Jl) 

o  a  k  o  v  ' 
o 


a-a 


(a  t ;T/2)  and  eUJ  2  .  f2 
o  '  a  k  o 

o 


is  negligible.  The  former  will  be  true  if 


a-a 


a  B  o 
o 


f  T  «T  and  the 


/  \  2 

latter  will  be  true  if  f2  T  tt  «  £  these  two  conditions  and 


a23 

o 


the  other  above  will  be  satisfied  if  (20)  of  the  main  text  is  satisfied. 
Comparison  of  the  two  results  above  yields 


xFM(at)  e 


a-a  f2 

o  «  0  0 

a-a 

a  k  XrM(a  t  +  — r^-  f  ) 

FM'  o  a  k  o' 
o 


(Eq.  E.3) 


On  the  other  hand  by  definition  (see  (lb)) 


(Eq.  E.4) 


a-1 


=  .l  XFMLt>T‘  a  V  xFMlat)  dt 


Substitution  of  Eq.  E.3  into  Eq.  E.4  yields 
a-a  f2 

D  “S-rrj— — —  -r-  a>  ,  a-a 

Y  (at, 1/a)  e  ao  /  xFM(t-T-  ^  tT)  xFM[a  (t+ — -  f  )]dt 

-oo  rn  a  I  rM  o  2.  o 

aoK 

a-a  f2 

-2nj  - r—  oo  a-a  , 

e  0(3  f  XFM^t-T"  ~7T  V  ~ a T  V  XFM^aot^  dt 


a2k 

0 


after  a  change  of  dummy  variable.  This  result  in  turn  yields  (19)  on 

comparison  with  the  equation  for  YrM(ant,  1/a  )  found  by  setting 

tT  =  0  and  a  =  a  in  Eq.  E.4.  0  0 

I  0  n 
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APPENDIX  F 


DERIVATION  OF  THE  BISTATIC  DOPPLER  SPREADING  FUNCTION  OF 
A  MOVING,  TURNING  LINE  TARGET 


The  derivation  takes  place  in  three  parts:  First  the  spreading  functions 
to  and  from  a  point  target  with  constant  range  rate  are  determined.  Next 
using  these  results  and  trigonometry  the  spreading  function  of  a  point  with 
unit  reflectivity  on  a  moving,  turning  line  target  is  found.  Finally  these 
points  are  summed  weighted  by  appropriate  reflectivity  to  obtain  the  overall 
spreading  function. 


F.l  THE  SPREADING  FUNCTION  TO  AND  FROM  A  POINT  WITH  CONSTANT  RANGE  RATE 


Consider  a  point  whose  range  p  is  given  by 

P ( t )  =  Py  +  py  t, 

then  for  transmission  to  the  point 

Pt  Pt 

y0U)  ■  yt<t  -  —  -  —  t)  . 

Consider  also  the  spreading  function 
S?(i,a;  ii,ai)  -  6(t-tj  )5(ot-ai ) 


From  (2b) 

y0(t)  =  yjtt-cuxi  +  (aj-1 )t] 


(Eq.  F.l) 


Comparison  of  the  two  equations  relating  yQ  and  yj  shows  that 


Sfp (^  >0^ iPf  >Pf )  ~  ^  (^"^yp)^  ( CX— OCy p )  , 


(Eq.  F.2) 


where 


aTP  ‘  1_pT/c>  tTP  ~  c 

c-py 

is  the  Doppler  Spreading  function  of  transmission  to  a  moving  point  target. 
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Clearly  for  reception  from  the  point  (with  pT  and  pT  replaced  by 
PR  and  pR) 


PR  +PRt' 

y0(t'+  -V-  ) 


=  yI(t')  . 


Pp^Ppi- 

Replacing  t'  by  t  =  t 1  +  — - -  yields 

Pp  i  Pp  Pp 

y0(t)  •  yj{t  -  (l  +  -f)  1  -f'K1  +  r]  -  1]  t  . 

Comparison  of  this  result  with  F.l  shows  that  the  spreading  function 
is 


SpR ( t »PR *PR )  ( T_TpR) 6 (ct_otpR )  ,  (Eg*  F»3) 

where 

“pr  4  <’  +  - 1  -  r  ■  tpr  4 

F.2  THE  SPREADING  FUNCTION  TO  A  POINT  ON  ft  MOVING,  TURNING  TARGET 


Application  of  the  law  of  cosines  to  Fig.  9  yields 


p|  =  r|  +  a2  -  2  r-j-  a  cos  a  , 

PR  =  rR  +  °2  '  2  rR  a  cos(a’Y)  • 

But  a  «  rj  and  rR  (on  the  order  of  102m  compared  with  on  the  order 

of  lC^m  to  105m);  therefore 


pT  =  rT  /  1-2  2-  cos  a  +  —  rT  -  a  cos  a  , 

T  T  rT  r|  T 

Pp  %  Cp  ~  o  cos(a-y) . 


(Eg.  F.4) 


Differentiation  of  these  results  yields 
pj  =  r-j.  +o  sina  a, 


PR  =  rR 


+  o 


sin(a  -  y)  (a-y)  , 


(Eg.  F.5) 


48 


SAC  ANTCEN  SR-44 


where  from  Fig.  9 

rT  =  -(v  cosa  +  Vj  cos  , 
rR  =  -  [v  cos(a-y)  +  vR  cos  <|»R]  . 

Now  consider  the  doppler  spreading  function  Sp(x,a;cr)  of  a  point 
at  location  a  on  the  target  assumed  to  have  unity  reflectivity. 

If  the  point  is  assumed  to  have  ranges  py  +  Pyt  and  pR  +  pRt  from 

the  transmitter  and  receiver  respectively  Fig.  3  may  be  used  to  obtain 
Sp  ;  therefore  Aty,  which  is  the  time  of  arrival  of  the  signal  at 

the  target,  must  be  selected  sufficiently  small  in  magnitude  that  these 
assumptions  are  valid. 

Use  of  the  fundamental  relationship  given  by  Fig.  3  and  (5b)  on  Eqs.  F.2 
and  F.3  yields 


Sp(a,i ;a)  =  /  /  -ayp)6(a ' (t-t ' )-Typ)6(a' -apR)6(i '-TpR)da 'di ' 


-00  0 


6(apR  "“Tp)  ^PR^'W'  ttP^ 


tTP 

=  <S(a-apR  otyp}6(T-TpR-  - — )  =  6[a-ap(a)]6[i-Tp(a)] 

PR 


since  a6 (at )  =  6 ( t) . 


It  is  convenient  to  approximate  ap(o)  and  Tp(o), 
(F-3) 

«p(o)  -  (l-^)/(H^)  . 


From  (F.2)  and 


Since  Py 
ap  (o) 


and  pR  «  c  (two  orders  of  magnitude  at  least) 

%  1  -  ( Py  +  pR)/c  , 

=  1  -  {ry  +  rR  +  o  [sina  a  +  sin(a-y)  (a-y)]}  /c 
=  1  -  (i  +  3  <5  o)/c  , 


(Eq.  F.6) 
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from  (F.5)  and  (23). 

t  (a)  = 


Pp 

no 


<'-r> 


Also  from  (F.2)  and  (F.3) 

p_l+p_ R 
c  c 


Pf  ,  PR  _  PT+PR  “pj0)'1  PT 

a  (o)c  c  c  a  (a)  c 
P  P 

rT+rR~°^C0S  a+  C0S  ap(a)-l 

c  otp(a)  c 


a(T)-l  a 
ap(o)  c 


cosa 


from  (F.4).  But 


ap(a)-l 
cxp  (a) 


1 

otp(a) 


Since  a  is  small  compared  to  r^  the  higher  order  terms  in  this 

expansion  may  be  neglected  in  the  third  term  of  t(o)  to  yield  with  the 
aid  of  (23) 

ap(r)-l 

xp(a)  «  [(£-6a)/c]  -  (rT/c)  -  +  aSo/c  , 


F.3  THE  OVERALL  SPREADING  FUNCTION 


Now  let  the  reflectivity  of  each  point  o  be  f(o);  then  from  linearity 
the  spreading  function  of  the  target  is 

n  00 

SD(x,a)  =  /  f(a)  Sp(i,a;  o)  do 
-00  ' 

This  integration  must  be  performed  for  two  different  cases:  6  =  o  and 
8  j*  o.  For  8  =  o,  ap(o)  =  1  -  i/c  and  is  independent  of  a  from  (F.6); 
hence 
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S(x,a)  «  /  f{a)  6[t-(£  +  ,rT  — ^ — )/ c  +  ( 1 -a ) 6  a/c]  6(a-1+£/c)  do 
-00  c4 


=  Tttfsi  f  TT^ijs  6fT"  (£  +  rT  -~-)/c+a]6(ot-l+i/c)  da 

C-£ 


c  rT 

( 1  -a ) 5  f(  (T-aTs  ^  <S(a-l+£./c)  , 


which  agrees  with  (22)  since  -(a-l)/a  =  l/(c-l)  from  the  6  function. 
For  6  i  o 


00  ap(a)-l  5  r  Aa 

S  ( x  ,a)  =  /  f{a)  6{t-£/c+[  ( 1  -  a  )6ct  +  —  rg\ —  r.j.]/c}6(a-l+|+  da 

-00  p'  ' 


66  f^d-a-i/c)]  6{t-  c  +  g(l-a)(l-a-  c)  +  a  c  )  • 


If  we  solve  for  ^  (1-a-  £■)]  from  the  argument  of  the  6  function  and 
substitute  the  result  into  the  argument  of  f  we  get 


f(- 


-ct+£-  — -  r- 


f— ■)«(!- i+JO-a)(l-at|)4St3jI) 


Finally  from  the  property  that  | a | 6 (at )  =  6(t)  ,  (22)  results. 


51 


SACLANTCEN  SR-44 


APPENDIX  G 

DERIVATION  OF  THE  MATCHED  FILTER  OUTPUT  FOR 
THE  BISTATIC  ECHO  FROM  A  MOVING,  TURNING  TARGET  IN  A  NON-DI SPERS I VE  MEDIUM 


Substitution  of  (19)  and  (kC)  into  the  fundamental  relationship  of  Fig.  I 
and  (5b) 

a ' -a  f 2 

oo  oo  -2ttj ( - -)  t—  n  a'-a  f  ,  , 

*J(t»d  f  e  a°  YFM[ao(T'T,+  _7"ir+ai^tT)1/ao] 

d  -oo  o  a 


„  -ex  +£ - r—  rT  0  0  rT 

^f( - -a— k  S[a’-V  |  ♦  VV>]  da'  dT' 


•  TTi]6'f( 


(l-a)6 


where  aQ  =  1 -£/c .  Interchanging  the  order  of  integration,  making  a 

0  a  I  _  1  rT 

dummy  variable  change  from  i'  to  x'--  +  —  ,  and  reversing  the 

order  of  integration  again  yields 


<b  (  t  .  1 )  ' 


a'-a  f2 

77  >  r  d 

/  /  e  o  YC( 


" 1  t 


1  +  3':1  At 
c  “Hf 


a'-an  f 

+  x?~£ir  >•  1/a<>]  •  f 

o 

a 

St'  fj 

=0  r 


(r$r) 6  a’-v  fe da  dT' 


__k_  /  e 
( 1 -a ) 6  -» 


0  r'  -CT’  \ 

f{rror) 


*  yFM  {%  T"T''  c  +  ^ - ~o7 1  )AtT+  — — - -  »l/a  }  dx1 

FM  C  H-a)a2  k  J  0 
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But,  as  pointed  out  in  the  main  text,  Aty  must  be  small;  therefore  the 
approximation 


_L .  Bt  1 
%  "  (l-a)o§ 


may  be  used.  If  it  is  used  and  the  dummy  variable  x'  replaced  by  Ax' 
where 


A  =  1  - 


B 

(l-a)a2. 


(Aty  +  ^  )  , 


the  result  is 


♦?<T-  ,!  *  77TT.  1  YFN[“o(t-T'-F  +  fj  stT>’  ■ 

3  A( 1 -a ) 6  -®  o 


.  e-*J  _ill  fA  f  (_2lL)  dT  '  . 
A(l-a)a0  k  A(l-a)6 


It  is  convenient  rewrite  the  last  result  in  terms  of  a  delayed  convolution 

4>^(x,  1)  *  ij>(x  -  Td), 

00 

^(t)  =  /  g(x  -  x' )  h(x' )  dx'  , 

-CO 


where  xd  is  given  in  Eq.  26  (recall  aQ  =  1  -  l/c)  and 


A(l-a)6 


A(l-a}6 
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9(t)  =  y“h  C<*0(t  -  V1  t2-).  '/<■„] 

0 


a  +1  a  -1  f 

V*  -  -s-  r  >  m  ,  , 

e  0  yFM  (T>ao} 


from  Eq.  12.  Fourier  transforming  yields 


=  F[*(t)] 


B  rFM  (f  -  fi’  “o>  H(f) 


where  ^(f.a)  and  f^  are  defined  in  Eq.  26  and 


H(f)  -  T*  Ch (x ) ] 


a  -1  f 

-2ttj  f'  — 

B  =  e  0  “o  k 


But  if  Y^(f)  and  HM(f)  are  defined  by 


H»(f)  =  b  f  (f  -  r) 


H(f )  =  Hn  (f  -  f‘) 


then  clearly  r^f)  obeys  the  equation  given  for  it  in  Eq.  26; 
furthermore  inverse  Fourier  transformation  of  ^(f )  yields 

ex  -1  f 

2irj  f '  (t  -  — —  r2-)  M 
vjj  (  t  )  =  e  %  iMt) 


(t,1) 

3 


27TJ  fO  ^T'Td  -  lh  JT}  M,  a 
;  o  i|»  (t  -  xd) 


which  agrees  with  Eq.  25  if  is  as  given  in  Eq.  26. 
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M 

To  obtain  H  (f)  it  is  necessary  to  Fourier  transform  h (t ) . 
from  Eq.  26 


Note  that 


-»-<  .  .  p  OT  o 

r  7TT1T 

[e  aoc  k  f(-t)] 


=  p  f  f  +  —  - _ —  i  • 

nr  c  a  k  >' 


therefore 
H(f)  = 


rrA(l-a)6  x  ,  6_  8 _ o  -I 

-  T  c  a  k  J 

o 


f  f2 

F(6{[1  -  a  -  ^  (AtT 


from  the  definition  of  A.  Frequency  shifting  this  result  gives  the 
M 

equation  for  H  (f)  given  in  Eq.  26. 
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